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1 Introduction 



Symmetries are one of the most fundamental concepts for understanding the laws 
of nature leading to conserving quantities. Unexpected violations of symmetries 
indicate some dynamical mechanism beyond the current understanding of physics. 

Parity violation was discovered in 1957 [[[]] in nuclear j3 decays and pion and muon 
decays 0. In the charged current interaction of the standard electroweak theory, 
parity and charge conjugation symmetries are maximally violated due to the V — A 
structure ||. All the experimental results up to now are in full agreement with the 
theory. 

A surprising discovery of the CP violating Kl — > 7r + 7r~ decays 0] was made in 
1964. The neutral kaon system still remains to be the only place CP violation has 
been seen. The Standard Model with three Fermion families can accommodate all 
the observed CP violation phenomena through the complex quark mixing matrix, 
Cabibbo-Kobayashi-Maskawa (CKM) matrix ||. However, no real precision test 
has been made due to the large uncertainties in evaluating the effect of hadronic 
interactions. 

Interest in CP violation is not limited to elementary particle physics. It is one of 
the three necessary ingredients to generate observed excess of matter over antimatter 
in the universe . The amount of CP violation which can be generated by the Stan- 
dard Model appears to be insufficient for explaining the observed matter-antimatter 
asymmetry in the universe [7|], giving a strong motivation to search for new physics. 

For CP violation in some B meson decay channels, the Standard Model can make 
precise predictions with little influence from the strong interactions. Those channels 
can be used to test the predictions quantitatively to look for a sign of new physics. 
Also in the B meson system, CP violation is expected in many decay modes. The 
pattern of CP violation allows us to make a systematic qualitative comparison with 
the Standard Model predictions. Therefore, it is now widely accepted that the B- 
meson system provides in future an ideal place for testing the Standard Model for 
CP violation @. 

In this article, we first derive the formalism [10| describing the particle antipar- 
ticle system, with and without CP violation. Three different mechanisms which can 
generate CP violation are clearly classified, together with experimental observables 
which identify contributions from the different mechanisms. Then, CP violation in 
the neutral kaon system is analysed in this formalism. After a brief discussion on the 
Standard Model description for CP violation in the neutral kaon system, we proceed 
to the neutral B meson system. Following the discussion on some Standard Model 
predictions, some thoughts are made how the situation could change if there exists 
new physics contributing in the B meson system. 
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2 Description of Particle Antiparticle System 



2.1 Basic Formalism 

Let |P°) and |P°) be the states of a neutral pseudoscalar particle P°-meson and 
its antiparticle P°-meson at rest, respectively. They have definite flavour quantum 
numbers with opposite signs: F — +1 for P° and F — — 1 for P°. Both states are 
eigenstates of the strong and electromagnetic interaction Hamiltonian, i.e. 

(H st + H em ) |P°) = m |P°) and (H st + H cm ) |P°) = m |P°) 

where mo and m are the rest masses of P° and P°, respectively. The P° and P° states 
are related through CP transformations. For stationary states, the T transformation 
does not alter them, with the exception of an arbitrary phase. While CP is a unitary 
operation, T is an antiunitary operation. 
In summary, we obtain 

CP\P°) = e idcp \P°) and CP |P°) = e~ i9cp \P°) 

T |P°) = e^ T |P°) and T |P°) = e^ T |P°) (1) 

where the ^'s are arbitrary phases, and by assuming CPT |P°) = TCP |P°) it follows 
that 

Since T is antiunitary, it follows that 

Tc = c*T 



where c is any complex number. If we define 

T\a) = \a), T\0) = 0) 
antiunitary operation has to give 



On the other hand, 



hence 



(oi\P) = (a 



(a\(3) = (a\ (T-'TIP)) = (a\ (t- 1 
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We can then conclude 

(a\(T^}) = 

i.e. when the T operator changes the direction of the operation, it must be complex 
conjugated. 

If strong and electromagnetic interactions are invariant under the CPT transfor- 
mation, which is assumed throughout this paper, it follows that tuq = tUq. 

Now we switch on an interaction, V, and the P can decay into final states f 
with different flavours (|AF| = 1 process) and P° and P° can oscillate to each other 
(| AF\ = 2 process). Thus, a general state \ip{t)) which is a solution of the Schrodinger 
equation 

* J^(*)> = ( H st + H cm + V) \m) (2) 

can be written as 

|^(t))=a(t)|P ) + 6(t)|P°)+Ecf(t)|f) 

f 

where the sum is taken over all the possible final states f and a(t), b(t) and Cf(t) are 
time dependent functions; |a(t)| 2 , |&(t)| 2 and |cf(£)| 2 give the fractions of P°, P° and 
f at time t respectively. Since the weak interaction is much weaker than strong and 
electromagnetic interactions, perturbation theory can be applied in order to solve 
equation [2| Also with the help of the Wigner-Weisskopf method, which neglects the 
weak interactions between the final states ||, and we obtain 

'i($M$)-(«-4)($) 

where the 2x2 matrices M and T are often referred to as the mass and decay 
matrices. 

The elements of the mass matrix are given as 

M ii = mo S ii + mj) + Yv(mm (4, 

f V m — ±j f J 

where V stands for the principal part and the index i = 1(2) denotes P°(P°). Note 
that the sum is taken over all possible intermediate states common to P° and P° for 

The elements of the decay matrix are given by 

r ij = 2ir y £(i\V\f)(f\V\j)S{m -Ef) (5) 
f 

The sum is taken over only real final states common to P° and P° for i ^ j. 



(( a \ T ) \P) 
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If the Hamiltonians are not Hermitian, transition probabilities are not conserved 
in decays or oscillations, i.e. the number of initial states is not identical to the 
number of final states. This is also referred to as the break down of unitarity. We 
assume from now on that all the Hamiltonians are Hermitian, i.e. 

Kt)| 2 + |&(t)| 2 + £M 2 = i, 
f 

and also 

M- ■ = m* r ■ = r* 

Clearly |a(t)| 2 + |&(t)| 2 decreases as a function of time, hence A is not Hermitian. 

Since the CP operator changes a particle state into an antiparticle state, the 
following relation can be obtained if V is invariant under the CP transformation, i.e. 
(CP)- 1 V CP = V: 

CP : \A X2 \ = \A 2X \ , A n = A 22 . 

Since the T operator induces complex conjugation, which is equivalent to interchang- 
ing a bra-state and a ket-state, the following relation can be obtained if V is invariant 
under the T transformation: 



T: \A 12 \ = \A 21 \ . 

By combining the two, we obtain for the CPT invariant case: 

CPT : An = A 22 . 

For a rigorous proof, equations 0, f| and || are used. 
It follows that 

•if An ^ ^22, i.e. M n ^ M 22 or r u ^ F 22 : 
CPT and CP are violated 

•if \A X2 \ ^ \A 21 \, i.e. sin (ip r - cp M ) ^ : 
T and CP are violated . 



where ^>u = arg(^i2) and ifr = arg(r 12 ). Note that CP violation in the mass 
and decay matrices cannot be separated from CPT violation or T violation. 

While there is no fundamental reason to respect CP and T symmetries, it can 
be shown based on only few basic assumptions that no self consistent quantum field 
theory can be constructed that does not conserve CPT symmetry Therefore, we 



restrict our further discussion to the case where CPT symmetry is conserved: i.e. 

Mn = m 22 = m, r n = r 22 = r 
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thus 



An = A 22 = A . 
Differential equation R| can be reduced to 



d 2 ait) , dait) / . , ,,\ . , 

+ 2 1 A + [A l2 A 21 - A 2 ) a{t) = 



(7) 



for a(t), and a general solution is given by 

a{t) = C + e~ iX+ + C_e~ iA - 

where C± are arbitrary constants which can only be defined by the initial condition. 
For b(t), we obtain 



b(t) 



A 
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. dait) . , . 
i— — - A ait) 
dt v ; 



which can be used once a(t) becomes known. 
Insertion of a(t) into equation [7| leads to 

X 2 ± - 2AX ± + {A l2 A 21 - ti' 

from which the eigen-frequencies are obtained as 



A± = A ± A 12 A 2 i = m± - -r± 



by solving where 
and 



m± = 3ftA± = M ± 3? ( ^li2 4u 



,1/2 



21 



,1/2 



(9) 

For an initially pure P° state, we have a(t) = 1 and b(t) — at t — 0, i.e. 
C+ = C_ = 1/2, and the solution is given by 



|P°(t)> = a(t)\P°)+b(t)\P°) 

= f + (t)\P°) + Cf-(t)\P°) 



1 + ICP 



|P+)e 



— j A+ t 



|P_)e' 



»; a_ t 



(10) 
(11) 



where 



/ ± (t) = i( e -^*± e - <A -*) 
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and C is 

(12) 

The two states |P + ) and |P_) are the eigenstates of A± and are given by 

|P±)= , 1 (|P°) ±C|P°)) • (13) 



1 + ICI 2 

For an initially pure P° state, we have 



TO) = |/-(0|P°) + /+(0|P°) (14) 



c 



i + IC 



(|P + )e- iA +*- |P_)e- iA -*) . (15) 



2C 

While P have definite masses and decay widths (as seen from equations |TT] and 



IB] ), P° and P° do not and they oscillate to each other (see equations [H] and II ). 
2.2 CP Conserving Case 

If V remains invariant under the CP transformation, from equations [l], % and |5] it 
follows that 

M12 = M 2 ie~ i2ecp = M* 2 e~ i2ecp 



thus 

and 

thus 



argMi2 = -^cp + n7r, 

r 12 = r 2 ie~ i2ecp = r* 2 e- i29ci 



arg r 12 = -^p + n'-K, 

where n and n' are arbitrary integer numbers. 
For (, we have 

where n" is an arbitrary integer number. The two mass eigenstates |P + ) and |P_) 
become CP eigenstates 

cp|p±) = ±(-ir"|p±) . 

The mass and decay width eigenvalues, equations || and || become 

m± = M ±{-l) n+n "\M 12 \ 



6 



1) 




Re 



3) 



Im 



4) 




Re 



Figure 1: Relative phase relations for M 12 , -Ti2, and CP transformation phase 9qp 
when CP is conserved: 1) CP = +1 state is heavier and decays faster, 2) CP = +1 
state is heavier and decays slower, 3) CP = +1 state is lighter and decays faster, 
4) CP = +1 state is lighter and decays slower. 



By examining various combinations of n, n' and n" , we can show that the following 
four possibilities exist: 

1. n=even, n'=even: CP = +1 state is heavier and decays faster, 

2. n=even, n'=odd: CP = +1 state is heavier and decays slower, 

3. n=odd, n'=even: CP = +1 state is lighter and decays faster, 

4. n=odd, n'=odd: CP = +1 state is lighter and decays slower. 

Figure |I] illustrates the phase relations in a pictorial way. The choice of n" does not 
alter the conclusion and n" = can be adopted without any loss of generality. In 
this case, |P + ) is the CP = +1 state. 



and 



r± = r±(-i) 



r 12 
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2.3 CP Violating Case 

Let us consider the time dependent decay rate for the initial P° decaying into a CP 
eigenstate f, given by |(/|V|P°(i))| 2 , and that for the initial P° decaying into f, given 
by |(/|^|P°(t))| 2 : 

~A 2 

R { (t) ex \f + (t)\ 2 + \f-(t)\ 2 + ™ 

R { (t) oc 



!/+(*)!' 



\f-(t)\ 2 + w ® 



(16) 
(17) 



where the instantaneous decay amplitudes are denoted by A{ = (f|V|P°) etc. and 
equations [l^ and [L4| are used. 

Since Rf(t) and R{(t) describe the CP conjugated processes to each other, any 
difference between the two is an clear proof of CP violation. As seen from the first 
terms of equations |I6| and |TT| , CP violation is generated if \Af \ ^ |Af|. This is called 
CP violation in the decay amplitudes. 

From the second terms of Rf(t) and Rf(t), it can be seen that CP violation is 
generated if \(\ ^ 1 even if there is no CP violation in the decay amplitudes. From 
equations [11] and it is clear that the oscillation rate for P° — > P° is different from 
that for P° — > P° if \(\ 1, thus this is called CP violation in the oscillation. 



The third term can be expanded into 




r + (t)f-(t) 




r + (t)f-(t) 



for -Rf(t) and 




£(*)/-(*) +T^Kf p r + (t)f-(t) 



ICI 



for -Rf(t). If CP violation in P°-P° oscillation is absent, the first terms are identical. 
Even in that case, if 




CP violation is still present. Since the process involves the decays of P° (P°) from 
the initial P° (P°) and decays of the P° (P°) oscillated from the initial P° (P°) into a 
common final state, it is referred as CP violation due to the interplay between 
the decays and oscillations. 

If CP violation in P°-P° oscillation is small, i.e. (\(\ — l) 2 << 1, we can derive 

|sin(y9 r - <p M )\ « 1 
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from equation [L2|, where ipr = arg/ 1 i 2 and <pm = argMi 2 as already defined, 
introducing |Aj-/ M | << 1 as 

fr - <Pm = n it - A r/M 

where n is an integer number, the following two approximations are possible: 
a) <fr — arg r 12 base 



By 



C 



2\M 12 \ 


|A 2 | 




4|M 12 


2 + 


IA2 


2 



1 



M±(-l) n \M 12 \ 

r±\r 12 \ 



'-\) n+1 + i 



2\M 



12 



A 



12 



b) v^m = argM 12 base 



+ 



2|M 12 | 


IA2I 




4|M 12 


2 + 


IA2 


2 



-l) n + Z 



I A 



121 



2\M 12 \ 



-l <PM 



m± 



M± \M 12 \ 

r±(-iy\r 12 \ 



(19) 



3 Neutral Kaon System 
3.1 Adaptation of Formalism 

Now we adapt the above developed formalism to the neutral kaon system. As de- 
scribed later, observed CP violation in the K°-K° oscillation is very small. The two 
mass eigenstates are called Ks and Kl with corresponding masses and decay widths 
referred to as ms, mL, A and 71 respectively and they are known to be mg < 
and Is > A- Therefore, M\ 2 and J12 is almost antiparallel to each other, thus 
n = 1 in equation |18[ 

Since the kaon decay properties are experimentally well measured, enough infor- 
mation is available to calculate r 12 from the data, as described in Section |3~5| . We 
therefore adopt the ipr base given in the previous section. 

It follows that 

C = (1 - 2e)e- ivr (20) 
where the small parameter e is given by 

2\M 12 \ \ 



\M l2 \ 


A2I 


sin (<fr ~ 


- ¥>m) 


4|M 12 


\ 2 + 


A2 


2 



l- 



\r- 



12 
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and 



|Ks) 
|K L ) 



1 



1 + |e| 2 
1 

1 + lei 2 



|K°) + (1 -2e)e-^ r \K°) 
|K°) - (1 -2e)e~ ilpr \K ) 



From the measured lifetimes [[12 

1 



and 



i.e. 



r s = — = (0.8934 ± 0.0008) x 10~ 1U s 

r, s 



r L = — = (5.17 ± 0.04) x 10" 8 s 



Ar = r s - r L = (1.1174 ± 0.0010) x 10 10 s^ 1 

and the mass difference, 



Am = m L - m s = (0.5301 ± 0.0014) x 10 10 hs' 1 



we obtain, 



and 



0.24966 ± 0.00004 





M 12 \ 


\r 


12 




4|M 12 


2 




r 


12 


2 



2|M 



1.2 



\r- 



0.9488 ±0.0026 . 
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(21) 
(22) 



Since the lifetime of Kl is much longer than that of Ks, it is possible to produce 
a Kl beam. Therefore, many kaon experiments have been done using Kl beams. 



3.2 CP Violation in Oscillations 

The CPLEAR experiment observed CP violation in the K°-K° oscillation by measur- 
ing the difference in the oscillation rates between K° — > K° and K° — > K°. The initial 
neutral kaons were produced by pp annihilations: pp — > K°K~7r + and — > K°K + 7r~, 
where the initial flavour can be defined by the charge sign of the accompanying kaon. 
Semileptonic decays were used in order to determine the flavour at the moment of the 
decay. Since the K° contains an s-quark (and K° an s-quark), K° (K°) can decay only 
into e + 7r - z/ (e~7r + T7) instantaneously. Therefore, the initial K° (K°) can produce the 
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Figure 2: Measured rate asymmetry between the initial K° decaying into e + 7r~z/ and 
the initial K° decaying into e~ir + V as a function of the decay time in units of ts by 
the CPLEAR experiment. The solid line is obtained by fitting a constant value. 

final state e~7t + T7 (e + 7r~z/) only through the K° — > K° (K° — > K°) oscillation. From 
the two measured time dependent decay rates, R e -(t) and R e +(t), an asymmetry 

= R_ e+ (t)-R c -(t) 
TU R c+ (t) + R e -(t) 



is constructed as shown in Figure |2|. Using equations [10], [L4| and |2(J, it follows that 

i-ICI 4 



A T (t) 



4 3fc 



i + ICI 4 

and from the measured A T (t) = (6.6 ± 1.6) x 1CT 3 (K 

|(| = 0.9967 ±0.0008 ^ 1 

is obtained exhibiting a clear sign of CP violation and T violation in the K°-K° 
oscillation. 

The parameter \(\ can also be measured from the semileptonic branching fractions 
of Kl by the lepton sign asymmetry: using equations and ^(| we obtain [|1^ 



Si = 



B(Ki 



IT V) 



B{K X 



B(K L - 

i-ICI 2 



^7r-z/) + 5(K L 
29fc 



i + ICI 2 

(3.27 ±0.12) x 10~ 3 
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where I can be e or /i and B stands for a branching fraction. 
Using all the measurements, we obtain 



9fe = (1.64 ±0.06) x 10~ 3 

and 

arge = (43.50 ± 0.08)°. 



3.3 CP Violation due to Decays and Oscillations 

Since the two-pion final state is a CP eigenstate with CP = +1, K L decaying into 
7r + 7r~ is a CP violating decay. This was indeed the first observed sign of CP violation. 
A commonly used CP violation parameter 7/ + _ is defined as 



(7l+7l-\V\K L ) 
(7T+7T-|V|Ks) 



(23) 



where equations [H], ^ are used and A + _ and A + _ denote the K° and K° — > 7r + 7r~ 
decay amplitudes respectively. 

The parameter rj + _ can be measured from the time dependent decay rates for the 



initial K and K into n n . From equations [11] and [L5| the two rates are given by 



R+-(t) oc 2 e ~ rs * + l r /+-| 2e_rL * + 2\rj + -\e- pt cos(Amt - + _) 



and 



R+-(t) oc 



l+4^e 



-r s t 



2|?7 + _|e r 'cos(Amt 



where + _ is the phase of ?7 + _ and r is the Ks-Kl average decay width. The 
second term is CP violating Kl decays and the third term is due to the interference 
between the K$ decay and CP violating Kl decay amplitudes. Figure [3] shows [14 
the measured i? + _(£) and together with the CP asymmetry defined as 



.4, 



, . R. 
•(* = ^ 



-(t)-R+-(t) 



R+-(t) + R+-(t) 



where the interference term is well isolated. At around t = lOrg, the Ks decay rate is 
reduced to the level of the CP violating K L decay rate, thus the asymmetry becomes 
very large. 
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Neutrol-koon decoy time [t s ] Neutrol-koon decoy time [t s ] 



Figure 3: The time dependent rate distributions for the initial K° (solid circles) and 
K° (open circles) decaying into tc + tt~ as a function of the decay time in units of Ts 
obtained by the CPLEAR experiment. The rate asymmetry is also shown. 

This direct comparison between the two CP conjugated processes illustrates an- 
other straightforward demonstration of CP violation in the neutral kaon system. 
From the asymmetry, the value of r] + _ is measured to be fl4| 

|?7 + _| = (2.264 ±0.035) x 10~ 3 , = (43.19 ± 0.60)° 

which leads to _ 

9 (c^-) = -(3.099 ±0.048) x 10~ 3 

exhibiting that CP violation due to the interference between the decay and oscillation 
is present. 

3.4 CP Violation in Decays 

The two-pion final state can be in a total isospin state of / = or / = 2. The 1 = 1 
state is not allowed due to Bose statistics. Using the isospin decomposition, we can 
derive the K° and K° decay amplitudes to 7r + 7r~ to be 

A + - = Wf(27r(/ = 0)|V|K°) + 7|(2tt(/ = 2)|V|K°> 



and 



A + - = J~(27r(J = 0)|V|K°> + JI(2k{I = 2)\V\K°) 
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Using CPT symmetry and the S-matrix, the K° and K° decay amplitudes can be 
related and it follows that 

A, 

where and a 2 are the K° decay amplitudes into 2n{l = 0) and 2n(I = 2) states 
due to the short-range weak interactions and 80 and 82 are the n-n scattering phase 
shifts for the 1 = and 1 = 2 two-pion configuration at = itir respectively. It 
is important to note that the two-pion scattering is totally dominated by the elastic 
scattering at the energy scale of the kon mass. Similarly for the 7r°7r° final state, we 
have 

V <i V <3 

Zoo = -^^e i ( 5 "+^p-^) + ^ a * e Hfe+^-^; 

As seen from the amplitudes, £?(Kg — > 7r°7r°)/.B(Kg — > 7r + 7r~) would be 0.5 if 
a 2 = 0. Since the measured ratio is ~ 0.46 |14|], we can conclude that |a2/ao| << 1- 
It follows that 

= (1 _ 2e') e _i ( 2 ^+^T-e C p) ( 2 4) 

A+- 

where the parameter e' is given by 




1 

71 



«2 
fl 



sin(y? 2 - Vo )e^ /2+5 *~ 5o) (25) 



and v^o, 2 = arga , 2 . 

As seen from equation CP violation in the decay amplitude, |Af_| 7^ |A H | 

is present if die' 7^ 0. From equation ^5], this is possible only if 



sin(y9 2 — ipo) 7^ and sin(<5 2 — 8q) 7^ . 

i.e. both the weak and strong phases have to be different for the 1 = and / = 
2 decay amplitudes. More generally, there must be two processes leading to the 
identical final state and both the strong and the weak phases must be different 
between the two processes in order to generate CP violation in the decay amplitudes. 
It should be noted that from the measured tt-tt scattering phase shift values, we 



have 115 



arge' = (43 ± 6)° 
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Using equations ^0] and |24|, it follows that 



c 



^ _ 2 e _ 2 e ') e _i (v^ + 2 ^o+6t-6'cp) 
1 - 2(e + e) -t( Vr + 2ip + 6 T - 6 CP ) 



where the approximation is made assuming that the phase difference between /12 and 

AqAq is small, which will be justified later. From equation [231, rj^ can be derived 

to be 

?7+- = e + i((p r + 2ip + T - 9 CP ) + e . 
Similarly the CP violation parameter for the 7r°7r° decay channel, 7700, is given by 

7/00 = e + i{(pr + 2y? + 0t~ #cp) - 2e' . 

Thus, we expect CP violation parameters to be different between the ir + ir~ and 7r 7r 
decay modes if e' 7^ 0. It has been shown by four recent experiments, NA31 flIE 



E731 m, KTeV PJ and NA48 19 



V+- 



Voo 



1.0127 ±0.0028 



i.e. CP violation in the decay amplitude is present in the neutral kaon system. If we 
neglect (ipr + 2ip + #t — # CP ), it follows that 



Hi 



1 

6 



V+- 



3.5 Phase of Decay Matrix 

As seen from equation [|, evaluation of A2 involves the decay final states which are 
common to K° and K°, which are 2n(I = 0), 2n(I = 2), = 1), 3n(I = 2) and 
3ir(I = 3) states: 

r 12 « J2 A Mi) A Mi) + Yj A *Mi) A Mi) ■ 
1=0,2 1=1,2,3 

The contribution from the decay amplitude to the 2n(I = 2) state is suppressed 
by the AI = 1/2 rule and the small measured value of e'. The contribution from 
the three-pion decay amplitudes are suppressed by r L /T s and the measured upper 
limits for the CP violation parameter for the 7r + 7r _ 7r and 7r 7r 7r° final states. In 
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conclusion, the phase of A 2 is essentially given by the phase of the A amplitude, 
and it can be expressed as 



<pr « &rgA* A = -2(p - # T + 9 C p 

so that 

\<p r + 2cp + e T -e CP \< o(io- 5 ) . 

Thus \<fir + 2<fio + 9^ — 9 C p\ << |e|, justifying the approximations made before. 
3.6 The Standard Model Description 



In the framework of the Standard Model |2(| , the short range contribution to K°-K° 
oscillation is obtained from the box diagrams (Figure f|) to be 



2 tl 

r] 1 a c S(x c ) + 2rj 2 a c a t E(x c ,x t ) + r] 3 a t S(x t ) 



M box = _G F f K B K m K m w 



12vr 2 

where Gp is the Fermi constant, /k, -Bk and are the decay constant, B parameter 
and mass for the K-meson respectively and is the mass of the W-boson. The 
QCD correction factors are denoted by 771 = 1.38 ± 0.20, r] 2 = 0.57 ± 0.01 and 
773 = 0.47 ± 0.04 and S and E are known functions of the mass ratios, x\ = m^/m 2 ^ 
for top (i=t) and charm (i=c). Note that 

S(x c ) « 2.4 x 10~ 4 , S{x t ) « 2.6, E(x c , x c ) w 2.2 x 10~ 3 (26) 



for m c = 1.25 GeV/c, m t = 174 GeV/c and mw = 80 GeV/c |T2|. The parameters 
a c and <r t are the combination of the elements of the Cabibbo-Kobayashi-Maskawa 
quark mixing matrix (CKM- matrix), 



Vc 



CKM 



' Kid Kis Vub 
Kd V cs V c b 
V Kd K s Kb 



o"c = KsKd* and a t = KsKd*- We adopt the following approximation of the CKM 



matrix using the parameters introduced by Wolfenstein |pl 



l-A 2 /2 A AX 3 (p~t7]) 

Kckm~| -\-iA 2 \ 5 v l-A 2 /2 A\ 2 j (27) 

AX 3 (l-p-tfj) -AX 2 -iAX 4 7] 1 

where where p = p(l — A 2 /2) and fj = 77(1 — A 2 /2). The parameter A is known from the 
light hadron decays to be 0.221 ± 0.002. From the B-meson decays, | K c b| = 0.0402 ± 
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K K K K 

Figure 4: The box diagrams contributing to the K°-K° oscillations. 



0.0019 and \V nh /V ch \ = 0.090±0.025 are measured @, giving A = 0.823±0.042 and 
Vp 2 + V 2 = 0-41 ± 0.11. The 5-parameter takes in account the difference between 
(0|ifwK ± ) and (f\H w \K°) where (0| is the hadronic vacuum state and (/| is the 
common quark states between K° and K°. The theoretical evaluations for this value 
vary between 0.5 and 1. 

In addition to M^£ x , there are large contributions from long range interactions 
Mf 2 , which are difficult to evaluate. Therefore, theoretical predication for M 12 = 
MJS* + Mff cannot be given. The long range interaction involves only the light 
flavours and its contribution to M\i is real in the CKM phase convention; the imag- 
inary part of M12 is generated only by the box diagram. Therefore we can derive 

, N SMi 2 2QM 1 b ° x 
sm(v9 M ) - 



M l2 \ Am 



In the CKM phase convention, r 12 can be approximated as real. Therefore, it follows 
that 



3fc 



2 Am 



Although there are considerable uncertainties to evaluate numerically this expression, 
the currently allowed range of the Wolfenstein parameters, A, A, p and r\ gives a 
consistent value of 3?e with the experimentally measured value. 

Prediction of e' requires an accurate evaluation of the phase difference between 
ao and ci2- For the a Q amplitudes, the tree, the gluonic penguin and the electroweak 
penguin diagrams contribute. Only the tree and electroweak penguin diagrams make 
contributions to the a 2 decay amplitude. All the penguin diagrams are shown in 
Figure ^. Not only the short range interactions, but also the hadronic matrix elements 
with long range interactions have to be evaluated in the calculations. This makes 
the numerical determination of e' very difficult. Within the theoretical uncertainties, 
values of e' calculated with the currently allowed range of A, A, p and r] are consistent 
with the data. 
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Figure 5: Gluonic and electromagnetic penguins contributing to the K° — > 27r decays. 



3.7 CP Violation in Rare Decays 

Experimental detection of Kl — > tt°vV is clearly very challenging. The final state is 
a CP eigenstate with CP = +1. Therefore, observation of this decay is a sign of CP 
violation. In the Standard Model, the decay is generated by penguin diagrams or 
box diagrams as shown in Figure |6|. 

Since the final state consists with only one hadron, long range strong interactions 
do not play a role and the decay amplitudes can be denoted as 

(7T Z/]7|# W |K ) 



G"7rQ ijij 



(7r°z/77|# w |K ) 



J 0CP-0T 



Unlike for the K° — ► 2ir decays, 
that we could have a situation 



arg a n o vV could be very different from O , so 



sin 



+ 2, 



» 



|sin (2 0^77 - 20o ) 
lei . 



The K L decay amplitude then becomes 



(ti ui7\H w \K l ) 



a w o. 



2 i \ a^ vV \ sin(2 <p n o vV 



1 _ 2e)e' i{2 ^°--~ 2M 
- 20 o ) • 




Figure 6: The box and penguin diagrams generating K° — > 7r°i/z/ decays. 
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Using isospin symmetry, the hadronic matrix element of the K° — > -n^vv decay 
amplitude and that of the K + — > 7i + e + u decay amplitudes can be related as 

(7T°|if w |K ) = (7r°|# w |K + > . 

This allow us to express the branching fractions for Kl — > it^vv using the branching 
fractions for K + — > 7T°e + ^ as 

B(K L -> ir°uV) -- 



\(7T°UV\H W \K L )\ 2 

r L 



T + I V-as \ 2 2tt 2 sin W 
r + 2tt 2 sin W 



3 x 10 



-11 



where X is a known function and 0w is the weak mixing angle. Since the hadronic 
matrix element is taken from the data, the theoretical uncertainties in this determi- 
nation is very small. Also the imaginary part of the amplitude is dominated by the 
short range interactions which can be reliably calculated. Therefore, the theoretical 
prediction can be considered to be clean. 

It is interesting to note that the CP violation parameter 

(7T Z/!7|y|K L ) 

7hfiw (Tflvv\V\Ks) 

as defined in the 2n case has l^o^l >> |e|, although the both final states have 
CP = +1. 

The current experimental measurement for this branching fraction is < 5.9 x 10~ 7 
with 90% confidence by the KTeV experiment f22fl , which is still far from the expected 
number. However, there are several proposals to observe the decays in the near 
future. 

4 B-meson System 

4.1 The Standard Model Description 
4.1.1 Some Elements of The CKM Matrix 

Among the nine elements of the CKM matrix, five of them related to the third 
generation play important roles in the B meson system: Vtd, Kit,, Vt s , V c h and Vtb- 
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Figure 7: Three elements of the CKM matrix, V t d, Kb, and Vts and the definitions 
of 0i, 03 and 5<j)3. 

In the approximation given in equation the phases of the five elements are given 
by 

arg V td = -0i, arg Kb = -03, arg V ts = <50 3 + n, arg V ch = arg V th = 

where 

Tj Tj 

0! = tan -1 , 03 = tan -1 — , 503 = tan" 1 \ 2 7]. 

l-p p 

Figure [7| shows the angles in p and rj planes. Note that 0i and 03 are often referred 
to as (3 and 7. Clearly <50 3 is very small, ~ 0.02. 

4.1.2 Oscillation Amplitude 

In the Standard Model, B-B oscillation is totally governed by the short range inter- 
actions, i.e. the box diagrams. Furthermore, only the top quark plays a role in the 
box diagram due to the large top quark mass (see equation ^6|) and the structure of 
the CKM matrix; 

Wb) =f ,_ aro^tb) . 1 1 

as seen from equation |2T| . 

Therefore, the off diagonal element of the mass matrix, M12 is given by |20 



Afi2 = ^B d 5(x t )(V; d Kb) for B d (28) 
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where /b cP B^ d and m-Q d are the decay constant, B-parameter and the mass of the 
B d meson. 

Similarly for the B s meson, we obtain 

Mi 2 = s —— 2 r)B s S{x t ){V ts V th ) for B s 

where /b s , Bb s and ms s are the decay constant, B-parameter and the mass of the B s 
meson. 

The phase of M 12 is then given by 



arg M12 



arg(K td K b ) 2 + n = 20! + 7T for B d 
arg(y t *y tb ) 2 + vr = -2503 + vr for B s . 



The parameter 7\2 can also be determined by taking the absorptive part of the 
box diagrams with charm and up quarks in the loops. The phase difference between 
M 12 and r 12 is given by 

argM 12 -argr 12 = 7r + -! (— Y rj x J (1 - P? + V 2 ' Bd (29) 



3 KrrihJ 



A 2 :B S 



i.e. sin(argMi 2 — argi~i 2 ) is small for Bd and very small for B s . Note that Mi 2 and 
I\ 2 are antiparallel. Therefore, the approximations for (, m± and i~± given on page | 
are valid with n = 1. Since we will rely on the Standard Model description of M 12 , 
and our experimental knowledge of the decay amplitudes is still limited, we adopt 
b) 4>m base. We refer the mass eigenstate with larger mass as Bh (B-heavy) and the 
other Bi (B-light) with their masses and decay width are given by: 

m h = M+\M 12 \, r h = r-\r 12 \ 

and 



mi = M- \M l2 \, /] = r + \r } 



12 



respectively, and B h (Bj) corresponds to P + (P_) defined in equation [L3. 
For both B<j and B s , we can now derive 



Ar 

Am 



r 



12 



M 



12 



3nm 



2 



2m 2 w S(x t ) 



5 x 1(T 3 for B d and B s (30) 



for m b = 4.25 GeV, m-w = 80 GeV and m t = 174 GeV, where Am and AT are 
defined as positive: 

Am = m b — mi, AT = f\ — i\ . 
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Using the measured values of Am = (0.464 ± 0.018) x 10 12 hs' 1 and the average 
lifetime r = 1/T = (1.54 ± 0.03) x 10~ 12 s for the Bj mesons, where t is the 
averaged decay width, it follows that 

Ar , 

-^^4x 10~ 3 for B d 

r 

and AT can be neglected in the decay time distribution for the Bd system. For the 
B s mesons, using the measured lifetime (1.54 ± 0.07) x 10 -12 s, it follows that 

-r- m 0.1 for B s . 

r 

The effect of AT is still not large, but can no longer be neglected in the decay time 
distributions. 

The small decay width differences of the Bd and B s systems do not allow to 
separate one mass-eigenstate from the other, which can be done for the kaon system 
by creating a Kl beam. Therefore, CP violation cannot be established by just 
observing the decays as in the case of Kl — > 2tt. We either have to compare the 
decay rates of the initial B° and initial B° states or measure the time dependent 
decay rates of at least one of the two cases, i.e. either initial B° or B°. 

Since Am = 2IM12L one can extract 



|^td| 2 = A 4 A 6 



Pf + f 



i.e. p and rj, from the measured B°-B° oscillation frequency Am^ using equation [28. 
However, theoretical uncertainties in calculating the decay constant and B-parameter 
are considerable and limit the accuracy on the extracted value of |Vtd| 2 . If the Bg-B° 
oscillation frequency Am s = 2|M^ 2 | is measured, |Vtd| 2 can be determined with much 
small uncertainty by using the ratio Am&/ Am s , due to better controlled theoretical 
errors in fB d / fs s an d B# d / ' B-q s . However, the frequency of the Bg-B° oscillation is 
expected to be > 1/A 2 = 20 times larger than that of the B°-B° oscillation and we 
still have to wait for sometime before it is measured. 



Since IM12/A2I << 1, C given by equation ^ can be further approximated as 



C 



1 1 ( A 2 

l — — <y 



2 \M 



'12 



where <fM — argM i2 as before. Seen from equation |29] and 0, the approximation 
|£| w 1 is accurate to 10~ 3 or better. 

Similar to the kaon system, CP violation (and T violation) in the oscillation can be 
measured from the time- dependent rate asymmetry between the initial B° decaying 
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into semileptonic final states with e + or R+(t) and the initial B° decaying into 
semileptonic final states with e~ or /i - , R_(t). The asymmetry is given by 



R+(t)-R4t) i-icr 



R+(t) + R.{t) l + |Ch 

which is a very small signal. 
From now on, we assume 



0(10-1) for B d and << 0(1(T 3 ) for B 



C 



-1<PM 



for both Bd and B s and Ar = for Bj. 

In summary, the two mass eigenstates are given by 



|Bh) 
|Bi) 



V2 
1 

72 



B) +e-^ M |B) 



B) 



-* <Pm 



B) 



and 



mi = mo + IM12I, rn\ = mo — IM12I, Am = — mi 
for Bd and B s . For the decay width, we have 

Ti = r h for B d 

ri = r + IA2I, r h = r - |r 12 |, Ar = /] - r h for b s 



4.1.3 Time Dependent Decay Rates 

Since AT is small in the B meson system, it is more convenient to derive the time 
dependent decay rate from the particle-antiparticle base rather than the mass eigen- 
state base. Using, equations [lC] and [14] the time dependent decay rates for the final 
state f can be derived as 



Rf(t) oc 
Rt(t) oc 



\AA 2 



\M 2 „-* t 



2|CP 



[I+(t)+I-(t)] 

[i+® - 



(32) 
(33) 



where P is the averaged decay time, t = (r+ + r_)/2, and Af is the instantaneous 
decay amplitude for the P° — > f decays. The two time dependent functions, I+(t) 
and I~(t), are given by 



I + (t) 
I_(t) 



Ar 



Ar 



(1 + \Lff) cosh — t + 23?L f sinh — t 
(1 — \L{\ 2 ) cos Ami + 2SLf sin Ami . 
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The parameter Lf is given by 



L f = C- 



where A { is the instantaneous decay amplitude for the P° — > f decays. 

The time dependent decay rate for the CP conjugated final states f CP are derived 
to be 



-Rfcp (t) oc 
i? f cp (t) oc 



|A fC p| 



-rt 



7f(t)+7 CT (0 



|A fCP | 2 |C| 2 



-rt 



/f(t)-/ CP (t) 



(34) 
(35) 



where A f cp is the instantaneous decay amplitude for the P° — > f CP decays. Two time 



dependent decay rates, I+ P (t) and Il r (t) are given by 



tCP, 



fCP 



n (*) 



1 + |L f CP | 2 ) cosh — - t + 23?L f UF sinh ==f- t 



Ar 



CP. 



J CP (t) = (1 - |Lp p | 2 )cosAmt + 2SL f CP sinAmt 



where the parameter, L\ , is given by 



^cp = 1 ^f CP 
f ~ C ^ f cp 

and A f cp is the instantaneous decay amplitude for the P° — > f CP decays. 

The decay rates Rf(t) and i? f cp (t) are CP conjugate to each other and so are Rf(t) 
and Rfcp(t). If there exists any difference between the CP conjugated processes, this 
is a clear sign of CP violation. 

The final state f can be classified into the following four different cases: 

I. Flavour specific final state (A { = A { cp = or A { cp — Af — 0) 

II. Flavour non specific final state 

Il-a. CP eigenstate (At = A { cp and Af = A { cp) 

Il-b. mixed CP eigenstate (Af = Afcp and Af = A { cp) 

II-c. CP non eigenstate 
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4.1.4 CP Violation: Clean Case 



The contribution to the B° decaying into J / if) Ks is dominated by the tree diagram 
with V^,K S . Although there exist some contribution from the penguin diagrams, the 
dominant penguin diagram contribution has the CKM phase V^^ts which is close to 
that of the tree diagram (Figure §). Thus, we can safely assume that there is no CP 
violation in the decay amplitude and the ratio of the B° and B° decay amplitudes is 
given only by the CKM part. By noting that CP(J/"0Kg) = —1 we obtain 



A(B ^ J/^K S ) 
A(B°^ J/^K S ] 



IKbKsKr s Kd| 2 



Using the formulae developed in the previous section, the time dependent rates for 
the initial B° decaying into J/^Ks, Rj/^Kgit), and that for B° decaying into J/^K S , 
Rj/i>K 3 (t) are given by 

-Rj/</> k s (t) oc e' rt + 3 L 3 /^ Ks sin Am tj 



Rj/ip K s 

which allow to extract 



(t) oc e 



-rt 



SLj /v ,k s sin Amtj 



3 Lj/ipKg 



A(B ^ J/^Ks) 



-9 



( Kd v tb Kt, ^cs K* Kd ) 5 



A(B° -> J/<0K s ), 
With the Wolfenstein parameterization, it follows that 



|VT d KbKbKsK*sK 



ud 



>J/i/>K s 



— sin 20i 



The same argument holds for the B s — > J/ip <fr decays and from the time dependent 
decay rates 

- / Ar Ar 

Ri/i>4>{t) oc e~" ( cosh— t + 2^L J/V)(A sinh — t + SLj/^sin Ami 

— - / Ar Ar 

Rj/ip<t>{t) oc e~ r * fcosh — t + 23? Lj/^ sinh — t - 3, L J/V ,^ sin Am t 
one can extract 



3L 



3 



Cx 



sin 



Note that we assumed in the calculation above that CP(J/ip <fi) = +1, i.e. the 
J /if) 4> state is in the lowest orbital angular momentum state of / = 0. If there exists 
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Figure 8: Tree and penguin diagrams contributing to the B° — > J/^Kg and B° — > 
J/ip (p decays. 

the I = 1 state with CP(J/ip(f>) = — 1, the measured ^Lj/^ will be diluted and 
the fraction of the CP = — 1 state must be experimentally measured. If there is the 
same amount of CP = +1 state and CP = —1 state, ^sLj/^,^ will vanish. 

An even cleaner decay channel is B° — > D* T 7r ± . There is only one tree diagram, 
b — > c + W + followed by W + — > u + d, which contributes to the B° — > D*~7r + decays. 
The same final state can be produced from the B° decays with another tree diagram, 
b — > u + W~ followed by W~ — ► c + d (Figure |). Therefore, the time dependent 
rate for the initial B° decaying into D*~7r + is given by 



R B *-(t) oc e~ rt 
where 



(1- 




to 


(1 + 


Lb*-tt+ 


2 ) 



1 + ; — cos Am t + : — sin Am t 

(1 + (Ld.-^+I 2 ) 



A(B° - D-7T+) 

-^D*-7T+ = C >< 



A(B° -> D*-tt+) 

The weak phase of A(B° D*-tt+) is given by KbK*d and that of ^( B ° ^ D*-tt+) 
by V^Vud- The phase of L^*-^ is then derived to be 

arg Ld*-^ = arg Kb - arg M 12 + Lf S 
= -03 + 20i + ip s 

where ips is a possible strong phase difference between the b — > u + W~ and b — > 
c + W + tree diagrams. 

CP conjug ated decay amplitudes of A(B° -> D*-tt+) and A(B° -» D*-tt+), i.e. 
A(B — > D* + 7r~) and A(B° — > D* + 7r _ ) respectively, are obtained by taking the 
complex conjugate of the weak amplitudes while the strong phase remains unchanged. 
Thus for D* + 7r~ we obtain 



i?D*+(^) oc e 



-Ft 
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B 



B 



d 



Figure 9: Tree diagrams contributing for the B° — > D* 7r + and B° — > D* 7r + decays, 
where 



rCP 



1 A(B° -> D*+7T- 
C X A(B° -»• D*+ 



7T~ 



and the phase of Ljy*- n+ is given by 

arg L^»- n + = - arg Kb + arg M 12 + y? s 

= 03 - 201 + if S 

From the two time dependent decay rates, we can extract 3 — 
Note that 



rCP 

-^D*-7T+ 



KbKd 



KbKd 



X 2 Jp 2 r] 2 « 1 



i.e. the effect we have to measure is small. 

The CP conjugated time dependent decay rate distributions are given by 



and 



i?D*+(£) oc e" 



Rn*-(t) oc e 



-rt 



1 + 



1 - |L CP _ I 2 ) 

1 D * ' cosAmt 



:i + i ^ 



sin Am t 



-rt 



(1- 


£D*-7r+ 


2 ) 


(1 + 


-^D*-7T+ 


2 ) 



29ij£)*- 7r + . 

cos Am t — j — sin Am t 

(1 + ILd.-^+I 2 ) 



which can be used to obtain the same information. 

A similar method can be used for the B° — > DfK ± decays to extract 03 — 2503. 
The effect is larger since 



I T I ~ 

I^D7K+I ~ 



v; h v us 
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4.1.5 CP Violation: Not So Clean Case 



The penguin contribution to the — > 7r + 7r~ decay was originally thought to be 
small and the decay would be dominated by the b — > u + W tree diagram. However, 
the discovery of -B(Bd — > K ± tt t ) > _B(Bd — > 7r + 7r~) indicates that the contribution 
of the penguin diagrams to the B d — > 7r + 7r~ amplitude should be ~ 20% or more. 

Due to the penguin contribution, the phase of the B° — > 7r + 7r~ decay amplitude 
deviates from that of V* h . Furthermore, CP violation in the decay amplitude could 
be present. Evaluation of those effects involves calculating contributions from dif- 
ferent diagrams accurately. Strong interactions may play an important role as well. 
Therefore, this decay mode may not be ideal to make precise determinations of p 
and rj from CP violation. 



4.2 Case with New Physics 

Decay processes where only the tree diagrams contribute should be unaffected by the 
presence of physics beyond the Standard Model. Therefore, \V c b\ and \V n h\ obtained 
from the semileptonic decays of B mesons would not be affected by the new physics 
and A and p 2 + rf can be obtained even if physics beyond the Standard Model is 
present. 

New physics could generate B°-B° and Bg-B° oscillations by new particles gen- 
erating new box diagrams. They could also generate a tree level flavour changing 
neutral current contributing to the oscillation. Since these contributions are through 
"virtual" states, they contribute to M12 with little effect on /12, i.e. 



M 



12 — Mf 2 M + m^ p , r 12 — j 12 



r- 



SM 



where M^ 1 and rf^ 1 are due to the Standard Model and M^ p is the contribution 
from the new physics. The measured Am is given by 2|M 12 | and can no longer used 
to extract \V td \ 2 due to M$f. 
Since 



r 



12 



12 



pSM 
1 12 



Am 



remains small, CP violation in the oscillation remains small as seen from equation [31 
Therefore, 



c 



-1<PM 



is still valid. However, note that 



argM 12 ^argMff. 
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Decay amplitudes from the penguin diagrams can be affected by physics beyond 
the Standard Model since new particles can contribute virtually in the loop. There- 
fore, the modes such as Bd decaying into 7r + 7r~, Y^it^ may have some contribution 
from the new physics. 

Since the decays Bd — > J/-0Ks and B s — > J/V0 are tree dominated, they are 
little affected by new physics. Therefore we have 

A(W^ J/y;Ks) _ A(B S °^ J/-00) 
A(B° -> J/^ K s ) A(B° -> J/^ 0) 

with the phase convention due to the Wolfenstein parameterization and 

£j/Vk s , 3/4,4, = Te~ iVM ■ ~ for B d -> J/^K S and + for B s -> J/^0 

and studies of the time dependent decay rates give argMi 2 . 

The B d — > D*7r and B s — > D S K decays are generated by only the tree diagrams 
and not affected by new physics. Therefore we have 

axgL D *- w + = -03 - argMi2 + y? s 

and 

arg L D , +7r - = 03 + arg M 12 + y? s 

and studies of the time dependent decay rates provide argM 12 + 03- Similarly studies 
can be done for B s — > D S K. 

By combining the measurements of Bd — > J/ip Kg and — > D*7r or B s — > J/-0 and 
— > D S K, the angle 03 can be determined even with presence of physics beyond the 
Standard Model. By comparing the result from B d and that from B s , consistency of 
the method can be tested. Since the phase of V u ^ is given by 3 and its modulus is 
measured from the semileptonic decay, p and rj can be extracted. Once A, A, p and 
i] are known, Mf 2 M can be calculated and from the measured Am and argM 12 , the 
new physics contribution M^ p is obtained. This can be used to identify the nature 
of the new physics contributing to the oscillation. 

4.3 Experimental Prospects 

A possible experimental programme for the study of CP violation in the B meson 
system and search for physics beyond the Standard Model can be summarised in the 
following steps: 

f . Determination of \V c b\ and |VJ,b| from semileptonic (and some hadronic) decays. 
2. Measurement of Am for Bd and B s 
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3. Measurement of SLj/^Kc, 

4. Measurement of Lj/^^, L D . T7r ± and L D ^ K± 

The first step has been made by ARGUS and CLEO at T(AS) machines and the 
four LEP experiments. BABAR and BELLE at the new asymmetric T(4S) machines 
and CLEO will improve the precisions on those determinations. Future improvement 
of theory is also an important factor. Half of the second step, Am(B d ) was done 
by ARGUS, CLEO, UA1 at the SPS Collider, the four LEP experiments, SLD at 
SLC and CDF at the Tevatron. For Am(B s ), we may have to wait for the next data 
taking by CDF, DO and HERA-B. The third step will be made by BABAR, BELLE, 
CDF, DO and possibly HERA-B by the year 2005. 

After the second step, four parameters of the CKM matrix are all defined within 
the framework of the Standard Model, e.g. A, A, p and r\. The third step provides an 
additional information tan -1 77/ (1 — p) within the framework of the Standard Model 
and consistency of the CKM picture can now be tested. 

As demonstrated in the previous chapter, if physics beyond the Standard Model 
exists, the fourth step is needed to clearly establish the evidence of new physics and 
separate the effect due to the Standard Model and that from new physics. After the 
third step, only p 2 +rf will be known from \V n \,\ and the information on tan -1 77/ (1— p) 
is spoiled by new physics. Only after the fourth step, p and r\ can be determined, 
together with isolating the new physics contribution. 

For the last step, new generation of experiments with statistics much higher than 
10 10 B mesons are needed. The B s meson is an essential ingredient. After 2005, LHC 
will be the most powerful source of B mesons. Experiments must be equipped with 
a trigger efficient for hadronic decay modes to gain high statistics for the necessary 
final states. Particle identification is also crucial in order to reduce background. 
LHCb is a detector at the LHC optimised for CP violation studies with B mesons. 
The two general purpose LHC detectors, ATLAS and CMS can contribute only to 
a limited aspect of the fourth step. A proposed experiment at Tevatron, BTeV, can 
also make the last two steps. 

Clearly CP violation is expected in many other decay channels. For many of them, 
there are some theoretical problems for making accurate predictions. However, they 
can be used to make a systematic study which will provide a global picture whether 
CP violation can fit into the CKM picture. With all those experiments, we continue 
to improve our understanding of CP violation and hope to discover physics beyond 
the Standard Model. 
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